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Surma  ry 

An  infinite  thick-walled  pipe  contains,  and  is  surrounded  by,  an 
inviscid  fluid.  The  displacements  of  the  pipe’s  wall  satisfy  the  exact 
linear  equations  of  elasticity,  and  the  interior  and  exterior  fluids 
satisfy  the  scalar  Helmholtz  wave  equation.  The  differential  equations 
are  solved  by  means  of  Fourier  transforms  to  give  the  dispersion  relation 
connecting  frequency  with  axial  and  circumferential  wavenumbers.  The 
dispersion  relation  is  solved  numerically  to  give  examples  of  ’real’ 
axial  wavenumber  versus  frequency  plots  for  selected  circumferential 
harmonics.  - 
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INTRODUCTION 


r 


Interpretation  of  the  vibration  and  sound  radiation  of  fluid-filled 
pipes  is  facilitated  greatly  by  the  availability  of  axial  wavenumber 
versus  frequency  plots.  Fuller  and  Fahy  [1]  discuss  the  physical  signi¬ 
ficance  of  the  real  and  complex  branches  of  the  wavenumber-frequency 
plots  of  fluid-filled  pipes.  Walter  [2]  uses  plots  of  the  real  branches 
to  demonstrate  that  the  maximum  responses  of  an  air-filled  pipe  wall 
occur  at  frequencies  at  which  the  in-vacuo  pipe  and  rigid-walled  duct 
mode  wavenumber^-frequency  plots  cross  each  other.  James  [3]  shows  that 
the  sound  radiation  from  fluid-filled  pipes  is  better  understood  when  the 
real  branches  of  the  wavenumber^frequency  plots  are  available. 

Each  of  the  above  authors  uses  a  shell  theory  [  4)  to  describe  the 
pipe  wall  vibrations.  However,  much  of  the  pipework  in  industrial  use  is 
thick-walled,  so  it  is  necessary  to  compute  wavenurrber-frequency  plots 
based  on  exact  linear  elastic  theory  in  order  to  check  the  validity  of 
shell  theory.  Kumar  [5]  studies  the  real  and  complex  wavenumber  plots  of 
fluid-filled  pipes  using  the  exact  theory  for  axisyrrmetric  vibrations. 
Gazis  [9]  uses  exact  linear  theory  to  investigate  three-dimensional  free- 
wave  propagation  in  hollow  elastic  cylinders. 

Contained  in  this  report  is  the  mathematical  analysis  needed  to 
obtain  the  dispersion  relation  of  a  thick-walled  pipe  that  contains  and 
is  surrounded  by  a  fluid.  Exact  linear  elastic  theory  for  an  isotropic 
circular  pipe  undergoing  arbitrary  motion  is  used.  The  purpose  of  the 
work  is  twofold.  First,  it  provides  the  initial  stage  in  the  analysis  of 
wave-propagation  and  sound  radiation  by  layered  pipes,  corresponding  to 
previous  work  on  layered  media  [6]  .  Secondly,  it  is  a  computational  tool 
for  studying  free-wave  propagation  in  fluid-filled  thick  pipes  and  it  is 
particularly  valuable  for  making  comparisons  with  shell  theory. 

2.  PROBLEM  FORMULATION 


An  infinite  thick-walled  circular  pipe  contains  and  is  surrounded  by 
acoustic  fluids  of  possibly  different  density  and  sound  velocity. 

Figure  1  shows  the  geometry  of  interest.  The  time  harmonic  factor 
exp(-icot)  is  omitted  from  all  equations. 


It  is  assumed  that  the  displacements  (u^, u^.u^)  of  the  pipe  are 

governed  by  the  exact  linear  equations  of  elasticity.  It  is  convenient 
to  represent  these  displacements  as  Fourier  transforms 
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The  equations  of  elasticity  are  coupled  in  terms  of  the  variables 
ur#  and  u^,  so  it  is  desirable  to  express  these  displacements  in 

terms  of  unknowns  (F,G,H)  which  satisfy  wave-equations.  These  functions 
are  also  represented  by  Fourier  transforms 
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The  stresses  in  the  pipe’s  wall  have  the  Fourier  transforms 
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The  Fourier  transform  representations  facilitate  solution  of  the 
wave-equations  and  enable  relations  between  unknowns  to  be  presented  in 
'spectral'  form.  The  variable  a  is  the  axial  wavenumber,  and  n  is  the 
number  of  wavelengths  around  the  circumference. 


The  procedure  for  solution  is  as  follows.  First,  matrix  relations 
between  the  surface  spectral  stresses  and  displacements  are  obtained. 
Section  3.  Secondly,  the  spectral  pressures  due  to  the  interior  and 
exterior  fluids  are  expressed  in  terms  of  the  spectral  displacements  of 
the  inner  and  outer  surface  respectively.  Section  4,  Finally,  the  pipe- 
fluid  boundary  conditions  are  applied  and  the  conditions  necessary  for 
free-wave  propagation  are  established.  Section  5. 


3.  THE  ELASTIC  LAYER 


Figure  1  shows  a  section  through  a  pipe  whose  inner  boundary  is  r=a 
and  whose  outer  boundary  is  r*=b.  The  surfaces  are  subject  to  prescribed 
normal  and  tangential  spectral  stresses 


r  r  w  r  b  b  b  a  a  a  ,1 
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which  produce  spectral  surface  displacements 

re  ,  b  b  b  a  a  a  .T 
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A  relation  is  required  between  the  suiface  spectral  stresses  and 
displacements. 


The  linear  elastic  equations  of  motion  [  7  ] 

2 

(Xni)grad(div  u)  +pVu=P32u/3t2 


0.3) 


are  reduced  to  three  wave  equations 

V2F+k2F=0  ,  92G+k2G=0  ,  ^H+k^O  0.4) 

1  s  s 

by  means  of  the  substitutions 
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uz  =  3F/3z+32H/3r2+(  1/r)  3H/3r+(  1/r2)  32H/3<|>2  =  3F/3z+V2H-32H/3z2 

Substituting  the  Fourier  transform  representations  of  F,  G  and  H 
into  the  equations  0.4)  gives,  after  integrating  the  resulting  equations, 
the  spectral  equations 
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The  spectral  displacements  are  obtained  from  equations  0.5)  and 
0.6)  as 
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uz(n,a,r)  = 

+(a2-k2)J  (Y  r)A  +  (a2-k2)Y  (y  r)Ac 
snss  sns6 

The  stress-displacement  relations  necessary  for  subsequent  analysis 
are 
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which  can  be  represented  spectrally,  after  simplifying  by  the  use  of  the 
Bessel  function  differential  equation 
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The  boundary  conditions 
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yield  the  matrix  equation 

[  P(n,a)]  [  A(n,a)]  =  [x(n,a)] 
6*6  6*1  6*1 


(3.10) 


and  the  spectral  displacements  equations  (3.7)  evaluated  at  r=b  and  r=a 
give  the  matrix  equation 

[  R(n,a)] [  A(n,a)l  =  [  u(n,a)]  (5.11) 

6*6  6*1  6*1 

The  elements  of  the  matrices  [P(n,a)J  and  [R(n,a)l  are  given  in 
Appendix  A. 

Equations  (3.10)  and  (3.11)  may  be  used  to  eliminate  [A(n,o)l  to 
give  the  formula 

(P(n,a)][R(n.a)]_1  {u(n,a)l  =  [x(n,a)]  (3.12) 

6*6  6*6  6*1  6*1 
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which  is  of  particular  usefulness  when  analysing  layered  media  by  the 
finite  element  method  [  6 ]. 

4.  INTERIOR  AND  EXTERIOR  FLUIDS 

(a)  Interior  Fluid 

The  linear  acoustic  equation  l8l 


V2Pi  =  (1/c2)32p./3i? 


(4.1) 


is  solved  by  replacing  p^tr.^.z)  by  its  Fourier  transform 

representation  and  then  integrating  the  resulting  equation.  Its 
spectral  solution  is 

Pi(n,a,r)  -  +  B2Yn(yir)  (4.2) 

The  finiteness  of  p^(n,a,r)  at  the  origin  r=0  requires  B2  to  be 

set  to  zero.  The  relation  between  the  fluid  pressure  and  displacement, 
viz 


3pi/3r=piu2ur 


(4.3) 


evaluated  at  the  boundary  r=a  enables  the  pressure  to  be  expressed 
in  terms  of  the  boundary  displacement  u^fn.a.a)  as 

p.  (n,a,r)  =  p^u2ur(n,a,a)  ^(y^rVy^tyja)  (4.4) 


(b)  Exterior  Fluid 

The  solution  of  the  wave-equation  relevant  to  outgoing  waves  is 

(8) 

Pe(n,a,r)  =  BgHp^r)  (4.5) 

The  boundary  condition  equation  (4.3)  Evaluated  at  i^b  enables  the 
pressure  to  be  expressed  in  terms  of  the  displacement  ur(n,a,b)  as 

p  (n,a,r)  =  pj2u  (n,a,b)H  (y  r)/y  H'(y  b)  (4.6) 

0  e  r  n  e  b  n  b 


5.  PIPE-FLUID  BOUNDARY  CONDITIONS 

The  boundary  conditions  to  be  applied  at  the  interface  between  an 
elastic  solid  and  an  inviscid  fluid  are  the  continuity  conditions  of  dis¬ 
placement  and  stress  in  the  direction  normal  to  the  interface.  The 
equations  (4.4)  and  (4.6)  reflect  continuity  of  displacement  and  the 
continuity  of  normal  stress  in  the  absence  of  external  forces  on  the 
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pipe’s  surfaces  is  given  by  the  equation 


[Pin*)]  [A(n,o)l 
6*6  6*1 


l-p  (n,a ,b) ,0,0,-p. (n,a,a) ,0,0] 
e  6*1  1 


(5.1) 


which  may  be  rearranged  by  making  use  of  equations  (4.4),  (4.6)  ana 
(3.11)  as 


[  P(n,a) +Q(n,a)H  A(n,a)  ]  =  0 
6*6  6*1 


(5.2) 


The  elements  of  [Q(n,a)  ]  are  given  in  Appendix  A. 

The  system  of  homogeneous  linear  equations  (5.2)  has  a  non-trivial 
solution  only  if  the  determinant  of  the  coefficients  vanishes.  For  given 
values  of  n  and  u,  there  will  be  real  and  complex  values  of  a  at 
which  the  determinant  vanishes.  Plots  of  these  values  versus  frequency 
are  called  wavenumber-frequency  or  dispersion  plots.  The  real  branches 
alone  are  of  interest  herej  they  are  the  wavenumbers  at  which  free-waves 
propagate.  The  complex,  or  purely  imaginary  branches  describe  the 
evanescent  waves  whose  effect  decreases  exponentially  with  distance.  When 
an  exterior  fluid  is  present,  there  will  be  no  real  branches  when  a<k 
(except  for  pure  torsion)  due  to  energy  loss  in  the  form  of  acoustic 
radiation  [  8  ]. 

6.  IMITERICAL  RESULTS 


Fortran  programs  have  been  written  to  compute  and  plot  the  real 
branches  of  the  axial  wavenumber  versus  frequency  plots.  A  root  is 
found  simply  by  stepping  through  a  range  of  ct-values  until  a  sign  change 
occurs  in  the  determinant:  it  is  refined  to  a  selected  accuracy  by 
repeated  interval  halving.  Care  must  be  taken  with  the  Bessel  function 
computations  otherwise  spurious  roots  may  occur.  Chosen  SI  constants  are: 

Steel:  E  =  19.5E10  o  =  0.29  p  =  7700.0 

Water:  p  =  1000.0  c  =  1500.0 

Figure  2  shows  the  real  branches  of  a  water-filled  pipe  whose  inner 
and  outer  radii  are  0.2096m  and  0.2350m  respectively.  The  frequency 
range  (2.8  times  'ring'  frequency)  and  thickness  to  mean  radius  ratio 
(0.114)  are  sufficiently  small  enough  to  allow  shell  modes  alone  in  the 
pipe's  wall.  The  physical  interpretation  of  the  plots  is  discussed 
elsewhere  ll  1„  The  plots  are  compared  with  those  obtained  from  a  shell 
theory  f 3 J .  The  significant  differences  at  the  higher  frequencies  occur 
when  the  waves  are  close  to  fluid-type  waves j  they  are  due  to  the 
assumption,  in  [3],  that  the  fluid  radius  is  equal  tc  the  mean  radius  of 
the  shell. 

Figure  3  shows  wavenumber  frequency  plots  of  an  in-vacuo  pipe  whose 
inner  and  outer  radii  are  0.10m  and  0.20m  respectively.  The  branches 


11 


labelled  4-10  involve  'thickness  vibrations'  of  the  pipe's  wall.  The 
plots  for  circumferential  harmonics  1  and  2  are  consistent  with  those 
obtained  by  Gazis  [  9  ]  who  discusses  the  nature  of  the  wavenumber  families. 
The  plot  for  n=0  has  many  features  in  common  with  a  plot  obtained  for  a 
0.10m  thick  plate  that  has  zero  wavenumber  for  the  anti-plane  strain 
branch  -  unpublished  work  arising  from  [6].  The  curious  behaviour  of  the 
branch  labelled  6,  for  n=0  and  1,  has  not  been  adequately  explained. 

It  is  also  present  in  plate  wavenumber  plots. 
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APPENDIX  A 


The  Matrix  [P(n,a)l 


11 

2yY23(;'(Ylb)  - 

*k23 

1  n 

O^b) 

12  = 

2yY2Y' ' (Y.b)  - 
In  1 

Xk2Y 

1  n 

(Y  jb) 

_ 

(2Un/b2  H  Y  b3 ' 

(Y  b) 

"  3n(Y 

b)] 

'13 

s  n 

s 

n 

s 

*14  = 

(2yn/b2H  Y  bY' 
s  n 

(Y  b) 
s 

-  Y  Or 
n 

sb)1 

Pl5  =-2i^Y2sJ[;i(Ysb) 

p16  =-2iayY2sY'*CYsb) 

p21  =  (2nu/b2][Tn(Ylb)  -  YlbJn'(Ylb)  1 

p22  =  (2nu/b2)[Yn(Ylb)  -  YlbY;(>xb)l 

p23  =  C y /b2 ) [  2YsbJUYgb)  +  (Y|b2  -  2n2)JMYsb)  ] 

p24  =  (u/b2)l2YsbY^CYsb)  +  [Y|b2  -  2n2)Yn(Ygb) ] 

p25  =  (2iaiin/b2)[  bYg3^CYsb)  -  ^(^b)  ] 

p26  =  (2iayn/b2)[bYsY'(Ygb)  -  YnCYgb)  1 

P31  =  2i«wYlJ^Ylb) 

P32  =  2iayYlY;(Ylb) 
p33  =  (iany/b) Jn(Ysb) 
p34  =  (iany/b)Yn(Ygb) 

p35  ’  "V2"2  ' 

p36  ■  u,s(2^  -  k|IY^tT8b) 

The  remaining  rows,  4,  5  and  6,  are  obtained  by  setting  b=a  in 
rows  1,2  and  3. 
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r11  •  Vn'^1 


r13  ■  (n/b)Jn(Ystp) 


The  Matrix  f  R(n,eQ  ] 


r 12  =  YlYn(Ylb) 
r14  =  (n/b)Yn(Ysb5 


r,,-  =  -i“Y  J'(Y  b) 
Id  s  n  s 


r1K  =  -i“T  Y'(Y  b) 
lb  s  n  s 


r_,  =  -(n/b)3  (Y  b) 
Z\  n  1 


r97  =  -  (n/b)  Y  (Y  b) 
ZZ  n  1 


r_-  =  -Y  3 '  (Y  b) 
Z3  s  n  s 


r9„  =  "Y  Y'(Y  b) 
24  s  n  s 


r_c  =  (i<*n/b)3  (Y  b) 
25  ns 


=  (i“n/b)Y  (Y  b) 
2b  ns 


r31  =  iaJn(Ylh) 


r32  =  iaVYlb) 


r33  =  0 


r34  =  0 


r35  =  (a2-ks)Jn(Ysb) 


r,,  =  (o2-k2)Y  h  b) 
3b  s  n  s 


The  remaining  rows,  4,  5  and  6,  are  obtained  by  setting  b=a 
rows  1,  2  and  3, 

The  Matrix  [  Q(n,a)  ] 

The  elements  of  [Q(n,a)l  are  identically  zero,  excepting  the 
following 


where 


’ll 

=  C  r,, 
e  11 

Q12  = 

Cer12 

°13  = 

Cer  13 

’l4 

=  Cer14 

°15  = 

Cer15 

°16  = 

Cer16 

!41 

=  Cir41 

°42  = 

Cir42 

°43  = 

Cir43 

'44 

=  Cir44 

°45  = 

Cir45 

°46  = 

Cir46 

Ce-> 

e“2Hn(Yeb)/YeH 

^Yeb) 

Ci  =  P 

i(,23n(Yia)/y.3 

'(y^) 

in 
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